Magnetic Domains and Stripes in the Spin-Fermion Model for Cuprates 
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Monte Carlo simulations applied to the Spin-Fermion model for cuprates show the existence of 
antiferromagnetic spin domains and charge stripes upon doping. The stripes are partially filled, 
with a filling of approximately 1/2 hole per site, and they separate spin domains with a tt phase 
shift among them. The stripes observed run either along the x or y axes and they are separated by a 
large energy barrier. No special boundary conditions or external fields are needed to stabilize these 
structures at low temperatures. When magnetic incommensurate peaks are observed at momentum 
7r(l, 1 — 5) and symmetrical points, charge incommensurate peaks appear at (0, 25) and symmetrical 
points, as experimentally observed. The strong charge fiuctuations responsible for the formation of 
the stripes also induce a pseudogap in the density of states. 

PACS numbers; 74.20.Mn, 71.10.Fd, 74.25. Ha 



In recent years neutron scattering experiments have es- 
tablished that magnetic incommensurability is a property 
common to most of the high-Tc cuprates |1[ . In addition, 
there is mounting evidence supporting charge stripe for- 
mation in these compounds as well, pj These nontrivial 
spin and charge arrangements may be crucial to under- 
stand the unusual transport and superconducting behav- 
ior of the cuprates. Early Hartree-Fock studies of the 
Hubbard model already predicted stripe formation with 
insulating characteristics. |^ Phase-separation between 
hole-rich and hole-poor regions in the Cu02 planes, sup- 
plemented by long-range Coulomb interactions, has also 
been proposed to explain the existence of stripes. Q In 
addition, ground states with metallic stripes made out of 
d-wave hole pairs have been observed in the t-J model 
using special boundary conditions. |^ However, a more 
detailed theoretical understanding of these phenomena 
and resolution of current conflicting results have been 
extremely challenging especially since the t-J and Hub- 
bard models used for the cuprates are considerably diffi- 
cult to study. As an alternative to this more traditional 
approach here we present the first computational study 
of a simpler phenomenological model, the Spin-Fermion 
(SF) model, which has been previously analyzed mostly 
using mean field approximations with the main goal of 
understanding d-wave superconductivity. [^-|| Here the 
focus is instead shifted toward the magnetic and charge 
properties of the SF model, which have not been explored 
before using unbiased techniques. In carrying out such 
a study unexpected results were observed, notably the 
presence of spin incommensurability and metallic stripe 
formation at finite hole density, in excellent agreement 
with experiments. 

The SF model is constructed as an interacting system 
of electrons and spins, mimicking phenomenologically the 
coexistence of charge and spin degrees of freedom in the 
cuprates. |6|-||]. Its Hamiltonian is given by 



H = -t 



h.c.) + J ^ Si • Si + J' ^ Si • Sj, 



(ij) 



where cj„ creates an electron at site 



(1) 
with 



spin projection a, Si=^^^ cl^(JapCip is the spin of the 
mobile electron, the Pauli matrices are denoted by cr. Si 
is the localized spin at site i, (ij) denotes nearest-neighbor 
(NN) lattice sites, t is the NN-hopping amplitude for the 
electrons, J > is an antiferromagnetic (AF) coupling 
between the spins of the mobile and localized degrees 
of freedom (DOF), and J' > is a direct AF coupling 
between the localized spins. The density (n)=l — x of 
itinerant electrons is controlled by a chemical potential 
fi. Hereafter t = 1 will be used as the unit of energy. 
From previous phenomenological analysis the coupling 
J is expected to be larger than t, while the Heisenberg 
couphng J' is expected to be smaller. 

Here classical spins with j^ij = 1 will be used for the 
localized spins, as also assumed in previous literature, 
This will allow us to perform Monte Carlo (MC) simula- 
tions of model Eq.(l) without "sign problems", reaching 
by this procedure temperatures as low as T=0.01 at any 
density. This temperature is well below T=0.2, the low- 
est that can be stabilized with quantum MC away from 
half-filling in the standard Hubbard model. JM] The value 
of J will be fixed to 2, as suggested in Ref. M where com- 
parisons with experimental results were performed. The 
coupling J' among the classical spins will be set to 0.05. 
This value was selected by monitoring the magnetic sus- 
ceptibility and comparing its behavior to experimental 
results for the cuprates. The present study has been per- 
formed mostly on 8x8 lattices with periodic boundary 
conditions (PBC), but occasional runs were made also 
using open and antiperiodic BC as well as different lat- 
tice sizes. The numerical technique used here involves a 



1 



standard Metropolis algorithm for the classical spins and 
an exact diagonalization for the itinerant electrons. The 
details of the method have been described in Ref. | |Tl[ |. 

To study the magnetic properties of the system we 
measured the spin-spin correlation functions among the 
classical spins defined as u;{r) = ^J^ii^i ' Si+r)j where 
N is the number of sites. The Fourier transform of a;(r), 
i.e. the spin structure factor S{q), was also investigated. 
The momentum takes the values 27rn/L^, with n run- 
ning from to — 1, and denoting the number of 
sites along the 7=x or y direction. In our MC simula- 
tions long-range AF order has been observed at (n) = 1.0 
as expected. As the electron density is reduced from 
1, the S{TT,Tr) intensity decreases. One of the main re- 
sults observed in this effort occurs at finite hole den- 
sity where a remarkable spin incommensurability ap- 
pears at (n)Ri0.8 with the structure factor peak moving 
to 7r((l — S),l) = (37r/4,7r) and rotated points. This 
behavior is illustrated in Fig.l-a where 5(q) on 8x8 lat- 
tices at T=0.01 is shown for different values of (n) along 
the path (0,0) - (0,7r) - (tt, tt) - (0,0). Results along 
(0, 0) - (tt, 0) - (tt, tt) - (0, 0) are in the inset. 
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FIG. 1. (a) Structure factor S'(q) for the localized spins 
versus momentum for 3—2, J'=0.05 and T=0.01 on an 8x8 
lattice and several densities along (0, 0) — (0, tt) — (tt, tt) — (0, 0) 
(inset: results along (0, 0) — (tt, 0) — (n, tt) — (0, 0))); (b) same 
as (a) but at T=0.05; (c) spin and charge distribution for a 
typical MC snapshot at (n)=0.9 and the same parameters as 
in (a). The arrow lengths are proportional to the local spin 
Sz(i) and the radius of the circles is proportional to the local 
density n(i), according to the scale shown. PBC are used. 

Note that the spin correlations along the two paths 
of Fig.l-a are different because at low temperatures the 
symmetry under lattice 7r/2— rotations appears sponta- 
neously broken. Similar ground-state properties but ro- 
tated by 7r/2 have also been observed in independent runs 
depending on the initial conditions jlj] for the classical 
spins indicating that there are two energy minima sepa- 
rated by a large barrier. At T=0.05 (Fig.l-b) and higher 
temperatures the rotational symmetry is restored but the 



intensity of the incommensurate (IC) peaks is consider- 
ably reduced (e.g. at (n)=0.75 the T=0.01 IC peaks are 
four times higher than at T=0.05). This may explain 
the low intensity of the IC peaks observed in the T=0.2 
Hubbard model quantum MC simulations. ]To[ | 

Experimentally it was observed that S «0.25 corre- 
sponds to the saturation value reached at hole density 
xw0.12, which persists up to xwO.25. ||l^ Due to the fi- 
nite size of the lattices studied here, incommensurability 
for values of 5 smaller than 0.25 cannot be detected ap- 
parently preventing us from comparing directly with the 
low-doping experimental results. However, we have an- 
alyzed typical spin configurations (snapshots) emerging 
from our MC simulations close to (n)=1.0 and we ob- 
served the existence of large AF spin domains in most 
of them, an example of which is shown in Fig.l-c for 
(n)=0.9. The existence of these domains clearly sug- 
gests that tendencies towards magnetic incommensura- 
bility appear in the system at x<0.20 as well. 

The origin of the short-range incommensurate mag- 
netic order in the cuprates is still not clear. Magnetic 
order due to charge order has been proposed as a pos- 
sible explanation. To explore the possibility of charge 
ordering in the SF model, A^(q) was studied here de- 
fined as the Fourier transform of the charge correlations 
n(r) = H J^iii'^i - (n))(ni+r - (n))), where Ui is the num- 
ber operator at site i for the itinerant fermions. In the 
AF phase at (n) = 1.0 N{q) was observed to present a 
broad peak at q = (tt, tt) due to negative charge cor- 
relations (charge repulsion) at very short distances, in 
agreement with previous calculations Q. As the sys- 
tem is hole doped the behavior of iV(q) becomes more 
temperature dependent. At T=0.05 the peak in N{q) 
remains at (tt, tt), but at T=0.01 a sharp peak appears 
at small momenta for 0.7< (n) <0.9 (Figs.2-a,b) indi- 
cating the existence of extended charge structures. If 
charge and spin incommensurability were related, stripe 
studies [|j predict that the peak in iV(q) has to appear 
at (2(5, 0) and symmetrical points. Then, incommensu- 
rability associated to 5=0.125 that cannot be explicitly 
detected in S{q) due to the size of our clusters can nev- 
ertheless be observed in the charge channel. Indeed in 
Fig.2-a the peak in iV(q) for density 0.87 is located at 
q = (7r/4,0) compatible with S ~ 0.125. At low temper- 
ature the peak can be observed along the x direction but 
not along y (see inset of Fig.2-a) indicating that in the 
MC runs described here the charge domains are along 
the y direction, causing the spontaneous breaking of ro- 
tational symmetry described before. This clearly can be 
seen in the charge distribution MC snapshot at (n)«0.85 
shown in Fig.3-a. Note that the holes are located along 
the magnetic domain boundaries. The presence of stripes 
was explicitly verified also using 12x12 clusters. The re- 
sult is in excellent agreement with experiments on 
Nd-doped LSCO where vertical stripes are observed at 
the Sr concentration x «:! 1/8. 

In some stripe scenarios p6| charge ordering is ex- 
pected to occur at higher temperature than magnetic or- 
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der. However, in model Eq.(l) both charge and magnetic 
ordering appear to occur at similar temperatures. Con- 
sider, for example, in Fig.3-b a typical MC snapshot at 
(n)=0.75 and T=0.01. The holes are here aligned along 
two vertical columns. S'(q) and A^(q) for this particu- 
lar snapshot are very similar to the averages shown in 
Fig.l-a and 2-a, clearly indicating incommensurate be- 
havior in spin and charge. As the temperature is raised 
to T=0.05, stripes are no longer observed (Fig.3-c) but 
there are still hole-poor magnetic domains which pro- 
duce the (small) incommensurate peak shown in Fig.l-b 
at this density. An equally weak feature appears in A^(q) 
but it is more difficult to distinguish because the back- 
ground raises with increasing momentum, rather than 
being flat (Fig.2-b). This may be the reason why the 
spin incommensurability is easier to detect than charge 
inhomogeneities in models such as Hubbard or t-J, where 
low temperatures are difficult to reach. As it can be ob- 
served in the snapshot of Fig.3-c even at T=0.05 hole-rich 
patches and antiferromagnetically spin aligned hole-poor 
domains coexist, suggesting that incommensurate charge 
and spin order occur simultaneously due to their mutual 
interactions in the SF model. 




(1,1/2) (1,1) (1/2,1(2) (0.0) (0,0) (1/2,0) (1,0) (l,l(2) (1,1) (l/2,l/2) (0,0) 

q q 

FIG. 2. (a) Structure factor A'^(q) for the electrons versus 
momentum for J=2, J'=0.05 and T=0.01 on an 8x8 lattice at 
several densities along the path (0, 0)-(7r, 0)-(7r, 7r)-(0, 0) (in- 
set: same as main figure but along (0, 0)-(0, 7r)-(7r, 7r)-(0, 0)); 
(b) same as (a) but at T=0.05. PBC are used. 

Our results arc in qualitative agreement with the con- 
clusions of Ref. where it was argued that stripes 
can be stabilized at realistic values of J/t in the t-J 
model without the use of long-range Coulomb interac- 
tions. Moreover, in the SF model here we showed that 
charge stripes can appear spontaneously without the 
need of using external staggered magnetic fields or spe- 
cial BC to pin them, jl^ The SF model provides a clean 
and easy to study framework for the analysis of stripe 
formation in models of correlated electrons. The origin 
of the stripes in our study can be understood in part 
by analyzing the behavior of (n) vs fi shown in Fig.4-a, 
where it is observed that the density changes rapidly be- 
tween 0.5 and 1. At (n)«0.5 there is a plateau indicating 
that this density is particularly stable. For 0.5< (n) <1 



substantial charge fluctuations are to be expected due to 
the large value of d(n)/(iyLt, involving regions whith den- 
sity close to 1 (spin domains) and to 0.5 (hole stripes) 
which correspond to the two highly stable densities that 
appear in the system. In fact, calculating the charge den- 
sity along the stripes we found that 0.5< (n)stripc <0.65, 
which seems to indicate the existence of approximately 
one hole every two Cu ions. The charge density on the 
hole poor regions, on the other hand, has a value very 
close to 0.9. These densities inside and outside the stripe 
are in excellent agreement with experiments . The SF 
model improves on early Hartree-Fock calculations that 
predicted a stripe density close to zero. 
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FIG. 3. (a) Spin and charge distribution for a MC snap- 
shot on an 8x8 cluster for J=2, J'=0.05, T=0.01, and at 
{n)~ 0.85;(b) spin and charge distribution for a snapshot at 
{n)=0.75 and the same parameters as in (a). The notation is 
as in Fig.l-c; (c) same as (b) but for T=0.05. PBC are used. 

It is important to note that the SF model does not 
phase separate in spite of its large compressibility. Using 
different lattice sizes and BC we have verified that there is 
no discontinuity in (n) vs. /i, while the energy presents a 
nearly straight line behavior in the density range between 
0.5 and 1.0 (Fig.4-b). Using starting MC configurations 
in which all the holes were together at the center of the 
cluster it was observed that this arrangement decays into 
stripes, and comparing the energy of the phase separated 
and stripe configurations the latter was found to have a 
lower energy than the former. 

As remarked before, stripe configurations were ob- 
served with several BC including periodic. One subtlety 
encountered in the latter is that for densities where a sin- 
gle stripe is stabilized, the PBC prevented the occurrence 
of a TT-shift in the spin domains because it would induce 
spin frustration (Figs.3-a). However, if PBC are replaced 
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by open BC (OBC) the stripe still appears and in this 
case a 7r-shift is observed (Fig.4-c). When the number of 
stripes is even, as in Fig.3-b, the 7r-shift is spontaneously 
obtained independently of the BC. This result is also in 
excellent agreement with experiments. JTsf 

Our simulations can produce dynamical information 
directly in real-frequency without the need of carrying 
out (uncontrolled) analytic extrapolations from the imag- 
inary axis. This is particularly important to compare 
theoretical predictions with the results of recent pho- 
toemission experiments on optimally doped LSCO which 
showed the development of a pseudogap at fj, as the tem- 
perature decreases. pq | In Fig.4-d we present the density 
of states (DOS) for (Ey=0.75 at T=0.05, 0.02 and 0.01. A 
pseudogap at uj=^ clearly develops for decreasing tem- 
peratures. This is a consequence of the strong density 
fluctuations discussed above, and it is similar to the phe- 
nomenon recently observed in the context of manganites 
where a pseudogap develops due to the coexistence of 
hole-rich and hole-poor domains. p9[ | 
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FIG. 4. (a) Density (n) vs ^ for J=2, J'=0.05 and T=0.01 
on an 8x8 lattice with PBC; (b) energy versus density for the 
same parameters as in (a); (c) spin and charge distribution 
for a MC snapshot at (n) ~ 0.85 and the same parameters 
as in (a) using OBC. The notation is as in Fig. 3; (d) den- 
sity of states as a function of w — /i at (n)=0.75 for different 
temperatures. The remaining parameters are as in (a). 

Summarizing, the SF model has been studied using 
MC techniques without a-priori assumptions on their 
properties. Magnetic and charge incommensurability has 
been observed upon hole doping in this model. The in- 
commensurability is due to the formation of AF domains 
separated by metallic stripes of holes. A 7r-shift is ob- 
served between the domains, and the stripes are partially 
filled with an electronic density of about 0.5. Charge and 
spin incommensurability appear correlated, and stripe- 



like configurations are obtained independently of the BC 
used and without the long-range Coulomb repulsion. The 
effect arises from the strong charge fluctuations between 
densities 0.5 and 1.0, which in addition produces a clear 
pseudogap in the DOS. Phase separation has not been 
observed in this study. The complex behavior of the SF 
model reported here, not anticipated in previous analysis, 
suggests that this model can be as useful for theoretical 
studies of the cuprates as the Hubbard and t-J models, 
while computationally it is considerably simpler. 
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